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In even-dimensional Euclidean space for integer powers of the Laplacian greater than or equal 
to the dimension divided by two, a fundamental solution for the polyharmonic equation has 
logarithmic behavior. We give two approaches for developing a Fourier expansion of this loga- 
rithmic fundamental solution. The first approach is algebraic and relies upon the construction 
of two-parameter polynomials. We describe some of the properties of these polynomials, and 
use them to derive the Fourier expansion for a logarithmic fundamental solution of the poly- 
harmonic equation. The second approach depends on the computation of parameter deriva- 
tives of Fourier series for a power-law fundamental solution of the polyharmonic equation. The 
resulting Fourier series is given in terms of sums over associated Legendre functions of the first 
kind. We conclude by comparing the two approaches and giving the azimuthal Fourier series 
for a logarithmic fundamental solution of the polyharmonic equation in rotationally-invariant 
coordinate systems. 
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1. Introduction 

Solutions of the polyharmonic equation (powers of the Laplacian operator) are 
ubiquitous in many areas of computational, pure, applied mathematics, physics and 
engineering. We concern ourselves, in this paper, with a fundamental solution of 
the polyharmonic equation (Laplace, biharmonic, etc.), which by convolution yields 
a solution to the inhomogeneous polyharmonic equation. Solutions to inhomoge- 
neous polyharmonic equations are useful in many physical applications including 
those areas related to Poisson's equation such as Newtonian gravity, electrostatics, 
magnetostatics, quantum direct and exchange interactions (cf. §1 in Cohl &: Do- 
minici (2010) [J]), etc. Furthermore, applications of higher-powers of the Laplacian 
include such varied areas as minimal surfaces llSj] , Continuum Mechanics [9| , Mesh 
deformation 0] , Elasticity [l3] , Stokes Flow \^ Geometric Design [iS] , Cubature 
formulae mean value theorems (cf. Pizzetti's formula) |il4] . and Hartree-Fock 
calculations of nuclei 12 il l. 
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It is a well-known fact (see for instance Schwartz (1950) (^Ifi], p. 45), Gel'fand &: 
Shilov (1964) (0], p. 202) that a fundamental solution of the polyharmonic equa- 
tion on (i-dimensional Euclidean space is given by combinations of power-law and 
logarithmic functions of the global distance between two points. In a recent paper 
(Cohl & Dominici (2010) [4]), we derived a complex identity which determined 
the Fourier coefficients of a power-law fundamental solution of the polyharmonic 
equation. The Fourier coefficients were seen to be given in terms of associated 
Legendre functions. The present work is concerned with computing the Fourier co- 
efficients of a logarithmic fundamental solution of the polyharmonic equation. One 
obtains a logarithmic fundamental solution for the polyharmonic equation only in 
even-dimensional Euclidean space and only when the power of the Laplacian is 
greater than or equal to the dimension divided by two. The most familiar exam- 
ple of a logarithmic fundamental solution of the polyharmonic equation occurs in 
two-dimensions, for a single-power of the Laplacian, i.e., Laplace's equation. 

We present two different approaches for obtaining Fourier series of a logarithmic 
fundamental solution for the polyharmonic equation. The first approach is algebraic 
and involves the generation of a certain set of naturally arising two-index polyno- 
mials which we refer to as logarithmic polynomials. The second approach starts 
with the main result from Cohl & Dominici (2010) and determines the Fourier 
series expansion for a logarithmic fundamental solution of the polyharmonic equa- 
tion through parameter differentiation. Series expansions for fundamental solutions 
of linear partial differential equations such as the polyharmonic equation are ex- 
tremely useful in determining Dirichlet boundary values for solutions on interior 
domains (see for example Cohl & Tohline (1999) [11). 

This paper is organized as follows. In ^ we introduce the problem. In ^ we 
describe our algebraic approach to computing a Fourier series of a logarithmic fun- 
damental solution of the polyharmonic equation. In 21 we give our limit derivative 
approach for computing the Fourier series of a logarithmic fundamental solution 
of the polyharmonic equation. In ^ we give some comparisons between the two 
approaches. In ^ we use the results presented in the previous sections to obtain 
azimuthal Fourier expansions for a logarithmic fundamental solution of the poly- 
harmonic equation in rotationally-invariant coordinate systems which parametrize 
points in d-dimensional Euclidean space. In Appendix [A] we present some neces- 
sary formulae relating to differentiation of associated Legendre functions of the first 
kind with respect to the degree. In Appendix [B] we present some of the properties 
of the logarithmic polynomials. 

Throughout this paper we rely on the following definitions. The set of natural 
numbers is given by N := {1,2,3,...}, the set Nq := {0,1,2,...} = N U {0}, 
the set of integers is given by Z := {0, ±1, ±2, . . .}, and the set Q represents the 
rational numbers. For ai, 02, . . . G C, if i, j G Z and j < i then Yln=i = and 
an = 1, where C represents the complex numbers. The set R represents the 
real numbers. 



2. Fundamental solution of the polyharmonic equation and the 
non-logarithmic Fourier series 

If $ satisfies the polyharmonic equation given by 



(-A)'=$(x) = 



(1) 
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where x G R'^, A : CP(R'^) — > CP^^(R'^), for p > 2 is the Laplacian operator defined 
byA:=^ + ... + ^,A;GN and $ G C^^(R'^), then <I> is called polyharmonic. If 
the power k of the Laplacian equals two, then ([T|) is called the biharmonic equation 
and <I> is called biharmonic. The inhomogeneous polyharmonic equation is given by 



(-A)*^$(x) = /,(x), 



(2) 



where we take p to be an integrable function so that a solution to ([2]) exists. 
A fundamental solution for the polyharmonic equation on R'^ is a function : 
(R*^ X R'^) \ {(x, x) : X G R'^} — )• R which satisfies the equation 



(-A)^0f(x,xO = c^(x-x': 



(3) 



for some c G R, c 7^ 0, where 6 is the Dirac delta function and x' G R'^. When 
c = 1, we call a fundamental solution of the polyharmonic equation normalized, 
and denote it by Q'^ : (R'^ x R'') \ {(x,x) : x G R"'} R. The Euchdean inner 
product (•, •) : R*^ x R'^ — )• R defined by (x, x') := xix'^ + . . .-irx^x'^, induces a norm 
(the Euclidean norm) || • || : R'^ — )• [0, 00), on the finite-dimensional vector space R*^, 
given by ||x|| := -y/ (x, x). In the rest of this paper, we will use the gamma function 
r : C \ —No — )• C, which is a natural generalization of the factorial function (see 
for instance Chapter 5 in Olver et al. (2010) l^^)- A fundamental solution of the 
polyharmonic equation is given by the following theorem. 



Theorem 2.1 Let d,k e'N. Define 



{k-iy. (k-d/2)\ 22fc-i7rrf/2 
r(d/2- A;)||x-x'||2'=-rf 



(A; - 1)! 22fc7rrf/2 



(log||x-x'|| - h~d/2,d) 

if d even, k > d/2, 

otherwise. 



where fip^d G Q is defined as fip^d '■= \ \Hp + -ffrf/2+p-i ~ -f^d/2-1] j Hj being 
the jth harmonic number 



^ 1 



i=l 



then Q'^ is a normalized fundamental solution for (— A)'^ on Euclidean space R"^. 

Proof. See Cohl (2010) [3] and Boyling (1996) 

A separable rotationally-invariant coordinate system for the polyharmonic equa- 
tion ([1]) on R'^ is given by 



Xl 




• ,^d-i)cos 






. ,^d_i)sin 




= xsi^i, ■ 


• • ,Cd-l) 



(4) 



Xd = XdiCl, ■ ■ ■,^d-i) 
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which is described by an angle (j) £ H and {d — l)-curvihnear coordinates 
('^1) • • • 1 Cd-i)- A separable rotationally-invariant coordinate system transforms the 
polyharmonic equation into a set of d-uncoupled ordinary differential equations 
with separation constants m G Z and kj for 1 < j < d — 2. For a separable 
rotationally-invariant coordinate system, this uncoupling is accomplished, in gen- 
eral, by assuming a solution to ([1]) of the form 

d~l 

^{x) = e'^^ 7^(6, . . . , id-i) n ^^(^i' m, A;i, . . . , kd-2), 

i=l 

where the domains of the functions TZ and Ai, for 1 < i < d— 1, and the constants kj 
for 1 < i < c? — 1, depend on the specific rotationally-invariant coordinate system. 
A rotationally-invariant coordinate system parametrizes points on the {d — 1)- 
dimensional half-hyperplane given by i;^ = const, and i? > using curvilinear 
coordinates (^1, . . . ,^d-i)- (For a general description of the theory of separation 
of variables see Miller (1977) [12].) The Euclidean distance between two points 
X, x' £ R'^, expressed in a rotationally-invariant coordinate system, is given by 



||x-x II = y/2RR' [x - cos( 

where the toroidal parameter x > 1; is given by 

d 



1 1/2 



i=3 

2RR' 



l\2 



X--= ^^^^ , (5) 



where i?, i?' £ [0,oo) are defined in (jj]) for x,x' G R'^. The hypersurfaces given 
by X > 1 equals constant are independent of coordinate system and represent 
hyper-tori of revolution. 

From Theorem 12. II we see that, apart from multiplicative constants, the algebraic 
expression : (R"' x R"^) \ {(x, x) : x € R^} — t- R of an unnormalized fundamental 
solution for the polyharmonic equation in Euclidean space R'^ for d even, k > d/2, 
is given by 

If (x,x') := ||x - x'f fc-^ (log ||x - x'll - /3fc_rf/2,d) . (6) 
By expressing if in a rotationally-invariant coordinate system we obtain 

"1 



[f(x,x') = (2i?i?') 



log {2RR')-(3p^, 



[x - cos((/) - (j)')Y 



+^ {2RR'Y [x - cos{4> - ^')Ylog [x - cos(</< - </>')] , (7) 

where p = k — d/2 G Nq. For the polyharmonic equation in even-dimensional 
Euclidean space R"^ with 1 < A; < d/2 — 1, apart from multiplicative constants, the 
algebraic expression for an unnormalized fundamental solution of the polyharmonic 
equation {)f : (R*^ x R^) \ {(x, x) : x € R'^} — ;> R is given by 



flf(x,x') 
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By expressing in a rotationally-invariant coordinate system we obtain 

f)^(x,x') = {2RR')-'' [;^_cos(0 (8) 

where q = 2k — d. 

By examining ([7]) and ([SD, we see that for computation of Fourier expansions 
about the azimuthal separation angle {(j) — 4>') of if. and f)^, ah that is required is 
to compute the Fourier cosine series for the fohowing three functions fxi^x ■ R- — 
(0, oo) and g^^- : R — )• R defined as 

fxi^) '■= ix - cosily , 
9xW ■= (.X - cos tpy log (x - cos xp) , and 
'■= (x-cosV')"^, 

where p E Nq, g E N and x > 1 is a fixed parameter. 

The Fourier series of is given in Cohl & Dominici (2010) [4] (cf. (4.4) therein), 
namely 

(. - cosV^)^ = (.2 - If/' i:e.cos(n^) ^-^|^"|^-,"^' p; (7^31) ' (9) 

where the Neumann factor e„ = 2 — 6n,o commonly occurs in Fourier series, 5n,o is 
the Kronecker delta, and 

n 

{z)n ■.= Yl{z + i-l), 
i=l 

for z £ C and n E Nq, is the Pochhammer symbol (rising factorial). We have 
used Whipple's formula in ([9]) (see for instance, (8.2.7) in Abramowitz & Stegun 
(1972) jl[) to convert the associated Legendre function of the second kind Qu ■ 
(1,00) — >■ C appearing in \d\ to the associated Legendre function of the first kind 
P!f : (1, 00) — )• R. The associated Legendre function of the first kind can be defined 
using the Gauss hypergeometric function, namely (Magnus, Oberhettinger Sz Soni 
(1966) [ui, p. 153) 

The Gauss hypergeometric function 2 -Fi : CxCx(C\-No)x{2; E C : |z| < 1} — > C 
can be defined in terms of the following infinite series 

77/1, \ {a)n{b)n n 

n=0 

(see for instance Chapter 15 in Olver et al. (2010) 
The Fourier series of is given in Cohl & Dominici (2010) [?] (Whipple formula 
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(8.2.7) in Abramowitz & Stegun (1972) [1] and cf. (4.5) therein), namely 



where G N. Since the Fourier series of h-^ is computed in Cohl & Dominici 
(2010) [J], we understand how to compute Fourier expansions of l)f ([8]) in separable 
rotationally-invariant coordinate systems. In order to compute Fourier expansion 
of if ([7]) in separable rotationally-invariant coordinate systems, all that remains is 
to determine the Fourier series of g-^. This is the goal of the next two sections. 



3. Algebraic approach to the logarithmic Fourier series 

Since x > Ij one may make the substitution x = cosh 77 to evaluate the Fourier 
series of g^.. For instance, it is given in the form of (coshry — cosV')^log(coshr7 — 
cos'i/'), where p G Nq. For p = the result is well-known (see for instance Magnus, 
Oberhettinger & Soni (1966) [lH, p. 259) 



log(cosh 1] — cos ip) = 1] — log 2 — 2 cos(n-i/'), (11) 

n=l 

which as we will see, should be compared with (|lUp for (7 = 1, namely 

^ 00 

— J]e„cos(n^)e-"^. (12) 



cosh 77 — cos lb sinh 77 

n=0 

Note that for 77 > we may write and therefore 7/ as a function of cosh 77 
since sinh 7/ = \/ cosh^ — 1, = cosh 77 + \/ cosh^ — 1 , and therefore 77 = 
log (cosh 7/ + vcosl?7/ — 1 ) . Now examine the p = 1 case for g^. If we multiply 



both sides of (jlip by (cosh 77 — cosip) and take advantage of the formula 

3(71-0) COS t/; = - |cos[(7i + l)-0] + cos[(ri — l)-0] I , (13) 



cos 



then we have 



(cosh rj — cos log (cosh rj — cos ip) = {ij — log 2) cosh ij 
— (7/ — log 2) cos ip — 2 cosh 7/ COs(7l7/^) 



n=l 



e ^ r. N -, e 



+ ^ C0S[(71 + 1)7/^] + ^ C0S[(7Z - 1)7/;]. (14) 



n ^ — ' n 

n=l n=l 
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Collecting the contributions to the Fourier cosine series, we obtain 
(cosh rj — cos ip) log (cosh 77 — cos ^) = {1 + rj — log 2) cosh rj 



sinh r] + cos ip [ log 2 — 1 — r/ — — e '^'^ 



^ e cos nijj , . -IN 
+ 2. — ^—^ (cosh r] + n smh rj) . 



n=2 



n{n? — 1) 



If we compare (|15p with (jlOp for g = 1, namely 



(cosh 7y — cos V')^ sinh^ 77 



^ 00 

— y e„ cos (n^)e~"^ (cosh r/ + n sinh r/), 



n=0 



we notice that the factor (coshr/ + nsinhry) appears in both series. 
For p = 2 m. Qy-, we use (jl3p and similarly have 

(cosh 77 — cos V')^ log(cosh r/ — cos ip) 

= (77 — log 2) cosh^ 77 — 2(77 — log 2) cosh 77 cos ij) 



+ (77 — log 2) cos^ ^ — (2 cosh^ ^ + 1) cos(7i7/;) 

^ n 

n=l 

+ 2 cosh 77 cos[(r7 + l)-0] + 2 cosh 77 cos[(n — 1) 



71=1 



- - J] ^ C0S[(71 + 2)7^] -i^^ C0S[(71 - 2)7/;]. 

2 ■^^ — ' 77. ^ — ' 

n=l 



n=l 



1 V— > e""'' 



2 ^ 77 

n=l 



If we collect the contributions of the Fourier cosine series, we obtain 
(cosh 77 — cos ip)"^ log(cosh 77 — cos t/^) = (77 — log 2) ^cosh^ ^ 2 ) 



+ 2 cosh 7? e"'' - -e"^'' + 
4 



-2(77 — log 2) cosh 77 — ( 2 cosh^ V ~^ 2^ ^ ^ 



+ cosh 77 e 



-2r? _ ig-3r, 



COS 7p + - (77 — log 2) + 2 cosh 77 e 



2 cosh^ 77 + l)e~^'' + - cosh 77 e"^'' - -e" 
3 



477 



cos 27/; 



-4E 



e '^^ cos {nip) r/ 2 \ 1 2 1 2 n 

[(n — 1) sinh 77 + 3n sinh 77 cosh 77 + 3 cosh 77J 



ra=3 



7l(n^ — l)(77-2 — 4) 
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By comparing (I17p with (llOp for q = 2, namely 
1 1 



(cosh 7/ — cos V')^ 2sinh^r7 

oo 

e„ cos(n^)e~"'' [(n^ — 1) sinh^ + 3n sinh rj cosh r/ + 3 cosh^ ry] , (18) 



oo 

X 

n=0 



then we notice that the factor ((n^ — 1) sinh^ rj + 3n sinh r/ cosh + 3cosh'^ r/) ap- 
pears in both series. We wiU demonstrate in why the identification mentioned 
in (flBI) and (fTSl) occurs. 



This algebraic approach for determining the Fourier series of will now be 
generalized. By starting with (llip and repeatedly multiplying by factors of (cosh rj— 
cosi/j), we see that the general Fourier series of g-^ can be given in terms of a 
sequence of polynomials Rp : (1, oo) — )• R, with p S Nq and A; G Z, as 

(cosh r] — cos ipY log (cosh rj — cos V') = {r] ~ log 2) (cosh r/ — cos ip)^ 



p oo 

e 



nrj 

+ 2 ^ (-l)'=+ii?^(coshr/) ^ cos[(n + k)^P]. (19) 

k=—p n=l 



We will refer to Rp as logarithmic polynomials with argument x = coshr/ (in our 
notation p and k are both indices) (See Appendix |B] for a description of some of 
the properties of the logarithmic polynomials). 

The double sum in ()19p is simplified by making the replacement n + k ^ n. 
It then follows that the resulting double sum naturally breaks into two disjoint 
regions, one triangular 

A := {{k, n) : —p < k < p — 1, k + 1 < n < p}, 

with p{2p + 1) terms and the other infinite rectangular 

B := {{k, n) : —p < k < p, p + 1 < n < oo}. 

By rearranging the order of the k and n summations in (jl9p , we derive 

(coshry — cos ^)^ log (cosh ?7 — cosV') = {rj — log 2) (cosh — cosV')^ 

p p-i 
+ J]cos(nV)e-"''r4'"-Hcoshr/) + ^ cos(nV')e"^ r:^'-"-\coshr/) 

n=0 n=l 

E°° cos(nV')e~"''' (vx / , X /^^N 
12 \ ( 2 2T^«.p(coshr/), (20) 
n(n^ — 1 • • • — ) 



where x^^p^ ,'Rn,p '■ (1,00) R are defined as 



r^;/=(coshr?) :=2 



-l)'^'+ie'=^i?^(cosh??) 



71 — A: 

k=ki 
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and 

3f?n,p(cosh7/) := ^J^^^^^^, r~P'P(cosh??), {n > p + 1), 

respectively. We can also write the Fourier series directly in terms of the logarithmic 
polynomials Rp as follows 

(cosh r] — cos ipY log(cosh r] — cos tp) = (f] — log 2) (cosh r] — cos il^)^ 

P n-l (_l)fc+lefcr,^fc(cosh7?) 

+2j;cos(nV.)e-"^ J] ^ 

n=0 k=—p 

^ -J^^ (-l)'=+ie'=''i?^(coshr?) 

-2 y cos (71^)6"" y ^ ^ p - 

^ ^ ^' ^ n + k 

n=l k=—p 

.2 ± cos,„^).- t '-'''"flf-""'. (21) 

n=p+l k=—p 

Note that by using ([9|), then we can express (coshry — cos'^)^ as a Fourier series, 
namely 

(cosh 7? - cos il^y = sinhP rj V e„ cos(raV^) ~ "^' p"(coth ??), (22) 

^ (p + nj! 

where p G Nq. 
If we define the function '^n,p '■ (l, oo) — ^ such that 

^(coshry) if n = 0, 

-Dn,p(cosh 77) + £'„_p(coshry) if 1 < n < p — 1, 

e-P''t-P'P-i(coshr?) ifl<n = p, 

o nT 7^ 5R„,p(cosh 77) if n > p + 1, 

— V ) ■ ■ ■ (n'^ — i)n 



q3„,p(cosh?7) := < 



{•n? — p^) ■ ■ ■ {n^ — 1) 
where Dn^p, -E„^p : (1, 00) — >• R are defined as 



i:>„,p(coshr?) 



and 



En,p{cosh.vi) 



e^^rl^;""-^ (coshry) if p > 2, 
ifpG{0,l}, 



e""''t-P'"~Hcoshr/) if p > 1, 
if p = 0, 
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respectively, then we can write (I20p as 



(cosh r] — cos ipy log(cosh ry — cos ^) = {r] — log 2) (cosh r] — cos ip)^ 

oo 

+ ^ cos(n^)q3„,p(cosh ry). 



n=0 



In fact, if we use (I22p . then we can express g^. as 

oo 

(cosh 7] — COS Tpy log(cosh T] — COS if)) = COS (n-0)^n,p (cosh rf) 



n=0 



+ (r/ - log 2) sinhP r? V e„ cos(nV') ^ ^^"^^ ,,''^' p"(coth r]). 
^ (p + n)! 

Furthermore, if we define Qn,p '■ (1, oo) — )• R as 

n„,p(coshr?) := q3n,p(coshr?) + ^ni-pUp - n)] _ 2) ginh^ ??P;(coth ??), 

(p + n j ! ^ 

then we can express as 

00 

(coshry — cos V')^ log(cosh ry — cosip) = cos(n'^)Qn,p(coshT/). (23) 

n=0 

4. Limit derivative approach to the logarithmic Fourier series 

We now use a second approach to compute the Fourier series for a logarithmic 
fundamental solution of the polyharmonic equation ([6]). We would like to match 
our results to the computations in ^ which clearly demonstrate different behaviors 
for the two regimes, < n < p and n > p + 1. By applying the identity 

d 

(coshry — cos^/')^log(coshry — cos^) = lim — (cosh ry — cos V')'^"''^, (24) 



where p E Nq, to 



(coshry - cos V)^ = sinh^ ry ^ ^ cos(nV^) pn^^^^^ ^25) 

n=0 ~^ 



where G C\-N (cf. (3.11b) in Cohl & Dominici (2010) [1]), one can compute the 
Fourier cosine series of (cosh ry — cos tlj)P log(cosh ry — cos ip), provided availability of 
the necessary parameter derivatives. 
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Applying ([24]) to ([25]), we obtain 



11 



(cosh rj — cos ipy log(cosh rj — cos ip) 



d 

lim — — sinh''"''^ ri 



n=0 



+ sinh^' r/ ^ (-!)"■€„ cos{niJj) 



n=0 



1 



lim -— T 

u-fO dv {v +P+1) 



d 



P"(cothr/) 



Note that for p G Z and n S Nq, the associated Legendre function of the first 
kind Pp vanishes if n > p+1- This is easily seen using the Rodrigues-type formula 
(cf. (14.7.11) in Olver et al. (2010) [1^]) 



and the fact that Pp{z) is a polynomial in z of degree p. The derivatives are given 
as follows: 



d 

lim — sinh*^^^ rj = sinh^ rj log sinh r?, 



(26) 



d 



1 



lim ^ , , , ,, 

u^o du [u +p + l)n {p + n)\ 



[tp{p+l) - ijj{p+l+n)] 



(27) 



and 



(28) 



where ip : C\ — Nq — )• C is the digamma function defined in terms of the derivative 
of the gamma function 



dz 



T{z) =: i^{z)T{z) 



(see for instance (5.2.2) in Olver et al. (2010) [3]). The degree-derivative of the 
associated Legendre function of the first kind in (j28p is determined using ()A2p and 
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(IA3|) . By collecting terms and using (f26l) . (f27l) . and (f28]) . we obtain 



(cosh rj — cos TpY log (cosh — cos ip) = {rj — log 2) (cosh — cos if^Y 

X [2-0(2^ + 1) - V(p + 1 + - V'(p + 1 - "-)] Pp"(cothr/) 

p-i 



+ (-l)Vsinh^,V^Ii^^ 



n=0 



p-n-i (-l)fe(2n + 2/c + 1) 



' ^ ' L ■ (2n+fc)!(p-n)!j / .v, ^ 
^„,+.(cothr?) 



-j^ _^ fc!(p+n)! 



fc=0 

P 



{p — n — k){p + n + k + 1 



+2(-l)VsinhP??^ 



-l)*^ cos(n'(/') 



^1 (p-n)! 

)^(2fc + l) 
(iJ- A;)(p + A; + 1)' 



k=0 



+2(-l)P+^p!sinhP7? ^ cos(n^/;)(n - p - l)!Pp~"(cothr/). (29) 

n=p+l 



5. Comparison of the two approaches 

The limit derivative approach presented in 21 might be considered, of the two 
methods, preferred for computing the azimuthal Fourier series for a logarithmic 
fundamental solution of the polyharmonic equation. This is because it produces 
azimuthal Fourier coefficients in terms of the well-known special functions, associ- 
ated Legendre functions. On the other hand, the algebraic approach presented in 
^produces results in terms of the two-parameter logarithmic polynomials Rp{x). 
As far as the author is aware, these polynomials are previously unencountered 
in the literature. By comparison of the two approaches we see how the logarith- 
mic polynomials Rp{x) (potentially a new type of special function) are intimately 
related to the associated Legendre functions. In this section we make this com- 
parison concrete. We should also mention that the following comparison equations 
resolve to become quite complicated as p increases, and they have been checked 
for < p < 10 using Mathematica with the assistance of an algorithm generated 
using ()B2p and ()B3P from Appendix [Bl 



By equating the Fourier coefficients using the two approaches we can obtain 
summation formulae which are satisfied by the logarithmic polynomials. For n = 



February 9, 2012 1:27 Integral Transforms and Special Functions logfour'arxiv 



Integral Transforms and Special Functions 13 

and p > 1 we have 

P (-l)'=+ie-*^'?i?^(coshr?) 

k=i 

= sinh^r/ [V'(2p+ 1) - il){p + 1)] Pp(coth7?) 

for l<n<p — 1, p>2we derive 

n — k ^-^ n + k 

k=—p k=—p 

{ i-lT 



[p + n)\ 



, Jzi]_ V ^ [ (2n+fc)!(p-n)!j «^ 

(p + n)! 4^ (p-n- A:)(p + n + fe + l^ ^n+fcl^uuiA/; 



-j^ _^ fc!(p+n)! 



fc=0 



for n = p, p > 1 we obtain 



fe=— p 



^[V^(2p+l)-V'(l)]pP(coth7?) 

^ (-l)'=(2fe + l) , ' 



and for n > p + 1, p > 0, we see that 

^ (-l)'=+ie'=''i?^(cosh7?) 
^ n — k 

k=—p 



(-1)P+V(^-P- l)!e"''sinhPT/Pp-"(cothr/). (30) 



We now have closed- form expressions for the finite terms given by (j20p . (|23p 
and (|20p in terms of associated Legendre functions of the first kind. We also have 
a proof of the correspondence for the "ending" function '^n,q-i mentioned in ^ 
Through ([30]), the function K„,p (cf. ([H]), (HZD, ([HD, and ^) is directly 
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related to the associated Legendre function of the first kind, namely 



3?n,p(coshr/) = 2(-l)P+ip!(p + n)!e"^sintf r/Pp-"(cothr?). 



Therefore through (jlOp we have 



; — = r5 — — y^e„cos(n^)e "^5}?„ „_i (coshr?), 

(coshTy-cosV)" 2[(g- I)!]2sinh29-i7y^ " "''^'^ 



which demonstrates the correspondences which was mentioned near (jl2p . ()16p . and 
(HHD for (EOl) and (|29D. 



6. Fourier expansion for a logarithmic fundamental solution of the 
polyharmonic equation 



Now that we have computed the Fourier series for g^, namely (pO|) (cf. ([23]) ) and 
(|29p . we are in a position to compute the azimuthal Fourier series for a logarithmic 
fundamental solution of the polyharmonic equation ([6|). 
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For instance, using (f22]) and (f29]) . we have 



[^(x,x') = 1 {2RR'Y [log {RR') + r? - /3p,,] (x^ - If/' 



n=0 



X [2V'(2p + 1) - V(p + 1 + n) - V'(p + 1 - n)] P^" | 



1 

--(2i?i?'rp!(x'-ir/'E 



cos[n((^ - (j)')] 
{p + n)! 



n=0 

p-n-1 (-l)fc(2n + 2A; + 1) 
fc=0 



1 + 



k\{p+n)\ 
(2n+fc)!(p-n)! 



(p — n — /c) (;? + + /c + 1) 



pn 



n=l 



{p — n)\ 



n-l 



-l)^(2/c + l) X 



(p-/c)(p + A: + l) VVx^, 



X 



X E cos[n(0-0')](^-p-l)!^'p"" ^^-^ 

n=p+l \VX ^, 
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Alternatively, using ([7]), ([9]), and (f2T]) we obtain 

[^(x,x') = i {2RR'Y [log {RR') + 7? - {x' - If'^ 



X 2^ e„ cos[n(0 — J r: R^ 



n=0 



n=0 



^1 + 

fc=— p 



n — k 



p-i 



n=l 



-{IRR'Y cos[n(0 - <p')] (x + VF^^ 
^1 



k=—p 



n + k 



oo 



n=p+l 



p (-i)'=+i(x + yF^) iJ^Xx) 

fc=— p 



n — k 



Using these formulae, we can, for instance, obtain the axisymmetric component 
of a logarithmic fundamental solution of the polyharmonic equation, namely 



t^(x,x') =\{2RRy{x^-lYl^ 



n=0 2 



[log(i?i?') + 77 - /3p,rf + 2V'(2p + 1) - 2i,{p + 1)] pA yj— y 



or 



^^(x,x') 



X 



= J(2i?i?')^(x' - 1)^/' [log(i?i2') + r? - /3,,,] 

n=0 2 \ _ 1 



fc=l 



A; 
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The above expressions for the axisymmetric component of a fundamental solution 
of the polyharmonic equation is one type of expression sought after in Tsai, Chen 
& Hsu (2009) [ii]. 



Appendix A. Derivatives with respect to the degree of certain integer-order 
associated Legendre functions of the first kind 

The derivative with respect to its degree for the associated Legendre function of the 
first kind evaluated at the zero degree is given in §4.4.3 of Magnus, Oberhettinger 
& Soni (1966) [IT}] as 



d_ 
dv 



J u=0 



(Al) 



An important generalization of this formula has recently been derived in 
Szmytkowski (2011) [13] • The degree-derivative of the associated Legendre function 
of the first kind for p,m £ Nq and < m < p (cf. (5.12) in Szmytkowski (2011) 
[isl ]) is given by 



d 

—P"^(z) 



P-(z)log^ 



+ [2V'(2p + 1) - ^(p + 1) - V(p - m + 1)] P^{z) 



p—m—l 



k=0 



+ (-1) 



{2k + 2m + l) 



1 + 



k\{p+m)\ 



m— 1 



{p — m)\ 



{p — m — k){p + m + k + 

2k + I 

{jp-k){p + 



(fc+2m)!(p -m)!_, 



(A2) 



and for m > p + 1 (cf. (5.16) in Szmytkowski (2011) [181 ]) there is 



d_ 
dv 



(-l)P+"+^(p + n)!(n - p - l)!P-"(2). 



(A3) 



Some special cases of (|A2p include for m = 



d_ 

dv 



Pu{z) 



J u=p 



Pp{z) log ^ + 2 [V(2p + 1) - V(P + 1)] Pp{z) 



k=0 



{p-k){p + k + l)' 
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for m = p 



Howard S. Cohl 



Z + l 



PP{z) log —— + [2V^(2p + 1) - V(p + 1) + 7] PS{z) 



p—1 



k=0 



{p-k){p + k + l) 



where 7 = is Euler's constant 7 ~ 0.57721566490. Of course we also have 

for m = p = 



du 



Pu{z) 



log 



u=0 



Z+l 

2 ' 



which exactly matches (|A1 



Appendix B. The logarithmic polynomials 

The logarithmic polynomials are nonvanishing only for —p < k < p and by 
construction, they are satisfied by the following recurrence relation 



(Bl) 



Prom (jlip we have that R^{x) = 1. This gives us the starting point for the recursion. 
We conjecture that the derivative of the logarithmic polynomials R^ is given by 



j-Rf^{x)=pR^\{x). 

It is evident by construction that these polynomials are even in the index k, i.e., 

R';{x) = R^\x). 
Some of the first few logarithmic polynomials are given by 

R^ix) = 1, 

R',{x) = x, Rt\x)=^- 
= l + x^, RfHx) = X, Rf{x) = ^, 

Rlix) = \x + x\ Rf{x) = 1 + \x\ Rf{x) = \x, Rf{x) = \. 

We can find the generating function for the logarithmic polynomials as follows. 
Let 



00 00 



F{x,y,z) = Y^ Pp(^)y''' 

p=0 k=—oo 
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k 



be the generating function for the logarithmic polynomials Rp. If we define the 
function 



Sp{x,y)= £ 

A:=— oo 

then using the recurrence relation for (]Bip we can show 

Sp{x, y)= {^+\{y + Sp-i{x, y). 
Combining this result along with the fact that -^0(2;) = 1, we have 

Sp{x, y) = + \{y + ~yj^ ' 



so therefore the generating function for the logarithmic polynomials R^ is given by 



F{x,y,z) 



1 



l-z[x + \[y + ^ 



An algorithm for generating the logarithmic polynomials can be obtained by 
solving the following set of difference equations 



«o(p) = 




ai{p) = 


iai(p 


a2{p) = 


|«2(P 



(B2) 



anip) = ^an{p - 1) + xan-i{p - 1) + \an-2{p - 1) , 

subject to the boundary conditions 

ao(0) = 1 
ai(l) = j;ao(0) 

02(2) = j;ai(l) + ao(l) \ 
an{n) = xan-i{n - I) + an-2{n - I) ^ 



(B3) 



where R^ = ap_\i^\{p) is given along diagonals for a fixed p — \k\. For instance, 
one can obtain R^^{x) = for p > 0, R^^^ ^\x) = -^^x^ for p > 1 and 
= + E^x^ for p>2, etc. 
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